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Summary 

The  analysis  procedure  of  an  existing  computer  pro¬ 
gram  for  predicting  the  off-design  performance  of  radial- 
inflow  turbines  was  extensively  modified  to  improve  its 
accuracy  and  to  extend  its  applicability.  The  rotor  inlet 
slip  factor  correlation  was  replaced  by  one  that  includes 
rotor  blade  sweep.  At  the  rotor  exit  a  sector  analysis  was 
used  in  place  of  the  original  mean-line  analysis  to  model 
the  large  radial  variations  occurring  there.  Disk  friction, 
clearance,  and  vaneless  space  losses  were  added  to  the 
loss  model,  and  the  method  of  computing  trailing-edge 
loss  was  changed.  Finally,  the  entire  analysis  procedure 
was  recoded. 

The  resulting  computer  program,  called  RTOD  (Radial 
Turbine  Off-Design),  predicts  the  performance  of  a 
radial  turbine  (with  or  without  rotor  blade  sweep)  as  a 
function  of  pressure  ratio,  speed,  and  stator  setting.  The 
program  models  the  flow  with  stator  viscous  and  trailing- 
edge  losses;  a  vaneless  space  loss  between  the  stator  and 
the  rotor;  and  rotor  incidence,  viscous,  trailing-edge, 
clearance,  and  disk  friction  losses.  The  stator  and  rotor 
viscous  losses  each  represent  the  combined  effects  of 
profile,  endwall,  and  secondary  flow  losses.  The  stator 
inlet  and  exit  and  rotor  inlet  flows  are  modeled  by  a 
mean-line  analysis,  but  a  sector  analysis  is  used  at  the 
rotor  exit.  For  a  pivoting  stator  the  leakage  flow  through 
the  clearance  gap  is  accounted  for. 

This  report  describes  the  program  modeling  changes 
and,  for  completeness,  presents  all  analysis  equations. 
The  program  input  and  output,  including  error  messages, 
are  also  described  and  are  illustrated  with  an  example 
case. 

Comparing  calculated  results  with  experimental  data 
shows  that  the  revised  program  predicts  mass  flow  rate 
better  than  does  the  original  program.  Calculated  results 
are  also  used  to  illustrate  the  potential  improvement  in 
off-design  performance  offered  by  rotor  back-sweep  for 
high-work-factor  radial  turbines. 


Introduction 

The  radial  turbine  has  been  used  extensively  in 
automotive  turbochargers  and  aircraft  auxiliary  power 
units  and  has  recently  been  given  serious  consideration 
for  primary  powerplant  application  for  both  automobiles 
and  rotorcraft.  In  automotive  and  rotorcraft  applications 
the  turbine  operates  over  a  wide  range  of  power  settings 
and  requires  accurate  off-design  performance  prediction 
for  a  successful  design. 

A  computer  code  for  the  off-design  performance  anal¬ 
ysis  of  radial-inflow  turbines  is  described  in  reference  1, 
and  additions  to  this  program  for  variable-area  stators 
are  reported  in  reference  2.  Although  this  program  has 


been  used  extensively,  it  still  requires  modifications  to 
extend  its  accuracy  and  applicability.  Current  interest  in 
highly  loaded  radial  turbines  has  resulted  in  designs 
having  rotor  back-sweep,  a  feature  that  this  program  can¬ 
not  accommodate.  Also,  the  analysis  model,  as  shown  in 
reference  1,  predicts  too  low  a  value  for  mass  flow  in  the 
region  of  rotor  choke.  In  addition,  the  loss  model  does 
not  include  all  of  the  losses  present  in  the  design  program 
of  reference  3,  and  the  trailing-edge  loss  does  not 
properly  account  for  the  effect  of  variable  stator  angle. 

To  overcome  these  deficiencies,  extensive  modifica¬ 
tions  were  made  to  the  analysis  model.  The  rotor  slip 
factor  equation  was  replaced  by  one  that  includes  rotor 
blade  sweep.  A  sector  analysis,  similar  to  that  of 
reference  3,  was  used  at  the  rotor  exit  in  place  of  the 
mean-line  calculation.  Disk  friction,  clearance,  and 
vaneless  space  losses  were  added  to  the  loss  model,  and 
the  stator  and  rotor  trailing-edge  loss  equations  were 
changed.  Finally,  the  entire  procedure  was  recoded. 

The  resulting  computer  program,  called  RTOD  (Radial 
Turbine  Off-Design),  which  is  described  in  this  report, 
computes  off-design  performance  of  a  radial  turbine  by 
modeling  the  flow  with  stator  viscous  and  trailing-edge 
losses;  with  a  vaneless  space  loss  between  the  stator  and 
the  rotor;  and  with  rotor  incidence,  viscous,  clearance, 
trailing-edge,  and  disk  friction  losses.  RTOD  can  analyze 
a  variable  stator,  including  the  leakage  flow  in  the 
clearance  gap.  Furthermore,  RTOD  determines  the  rotor 
exit  flow  variation  from  hub  to  tip  and  can  account  for 
the  effect  of  swept  rotor  blading.  The  user  must  specify 
the  stator  and  rotor  viscous  losses  (which  represent  the 
combined  effects  of  profile,  endwall,  and  secondary  flow 
loss)  to  match  a  reference  performance  point.  RTOD 
then  calculates  the  turbine  performance  over  a  specified 
operating  range. 

For  the  sake  of  completeness  this  report  describes  all  of 
the  modeling  assumptions,  but  with  emphasis  on  recent 
modifications,  comprised  in  the  computer  program 
RTOD.  The  program  input  and  output  are  described  and 
the  use  of  the  program  is  illustrated  by  example 
problems.  Results  obtained  from  this  program  and  from 
the  program  of  reference  1  are  compared  with  experi¬ 
mental  data.  RTOD  can  be  obtained  from  the  Computer 
Software  Management  Information  Center  (COSMIC), 
University  of  Georgia,  Athens,  Georgia  30602. 


Analysis  Modification 

This  section  describes  the  changes  that  were  made  to 
the  previous  analysis  method  and  loss  model  of 
references  1  and  2.  The  described  losses  are  applied  at  or 
between  the  flow  stations  shown  in  figure  1.  Station  0  is 
the  station  upstream  of  the  stator,  station  1  is  just  inside 
the  stator  trailing  edge,  and  station  2  is  immediately 


downstream  of  the  stator  trailing  edge.  Station  3  is 
upstream  of  the  rotor,  station  4  is  just  inside  the  rotor 
trailing  edge,  and  station  5  is  immediately  downstream  of 
the  rotor  trailing  edge. 

Rotor  Analysis 

The  rotor  analysis  was  extended  to  include  a  slip  factor 
correlation  that  can  account  for  swept  rotor  blading. 
Also,  at  the  rotor  exit  a  sector  analysis  replaces  the 
previous  mean-line  analysis.  The  details  are  as  follows 
(all  symbols  defined  in  appendix  A): 

Slip  factor.  -  In  reference  1  the  optimum  rotor  inlet 
tangential  velocity  is  calculated  from  the  Stanitz  slip 
factor  correlation 


obtained  from  the  velocity  diagram  work  by  subtracting 
the  losses  due  to  disk  friction  and  tip  clearance.  This 
means  that  the  sector  velocity  diagrams  are  not  consistent 
with  the  shaft  output  work. 

Loss  Model 

In  reference  1  the  losses  through  the  stator  and  the 
rotor  are  each  modeled  by  a  single  overall  loss.  In  RTOD 
a  vaneless  space  loss  (between  stations  2  and  3)  was 
added,  as  well  as  a  rotor  disk  friction  and  tip  clearance 
loss.  Furthermore,  the  method  of  calculating  trailing- 
edge  pressure  drop  was  changed.  These  changes  are 
discussed  in  this  section. 

Vaneless  space  loss.  -  The  program  of  reference  1 
assumes  no  pressure  drop  from  the  stator  exit  to  the  rotor 
inlet  (between  stations  2  and  3).  In  RTOD  the  total 
pressure  drop  is  obtained  from 


The  Stanitz  correlation  is  limited  to  radial  blading  only.  _  4flpV2 

RTOD  uses  the  Wiesner  slip  factor  correlation  of  ^  2 Dhgc 
reference  4,  which  is  a  function  of  blading  (sweep)  angle. 

The  correlation,  when  adapted  to  a  turbine  geometry, 

takes  the  form  where  /  is  the  Fanning  friction  factor  for  parallel  walls 

and  /  is  a  length  based  on  an  average  flow  angle  between 


A 


[l  -Vcos  abl/(^bl)0-7i[l  -  [(^■av/Q-<lim)/(l  ~  ^lim)|3]  \ 
1  -  tan  ab!/tan  a3  / 


For  radial-bladed  rotors  (no  blade  sweep)  the  difference 
between  the  two  correlations  is  quite  small.  For  the 
example  1  turbine  the  Wiesner  correlation  predicts  an 
optimum  inlet  tangential  velocity  2.7  percent  less  than 
that  predicted  by  the  Stanitz  correlation.  The  full 
incidence  loss  calculation  procedure  is  shown  by 
equations  (B71)  to  (B82). 

Exit  sectors.  -  The  program  of  reference  1  analyses  the 
rotor  exit  at  a  mean  radius  only.  In  RTOD  the  rotor  exit 
is  divided  into  a  specified  number  of  sectors  that  are 
coupled  by  simple  radial  equilibrium 

dp  _  p  Vj 
dr  rgc 

Temperature,  pressure,  and  velocity  are  calculated  for 
each  sector  (eqs.  (B94)  to  (B106)),  and  overall  values  of 
mass  flow  and  work  are  obtained  by  summing  the  indivi¬ 
dual  sector  values.  These  calculations  take  into  account 
the  viscous  loss,  since  a  viscous  rotor  loss  coefficient  is 
specified  for  each  sector.  The  disk  friction  and  tip 
clearance  losses,  however,  are  of  an  overall  nature  and  do 
not  reflect  hub-to-tip  variations.  The  shaft  output  work  is 
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stator  outlet  and  rotor  inlet.  The  calculation  procedure  is 
given  by  equations  (B45)  to  (B54). 

Rotor  disk  friction  loss.  -  The  program  of  reference  1 
does  not  calculate  a  separate  disk  friction  loss  but  lumps 
all  losses  (viscous,  disk  friction,  clearance,  and  trailing 
edge)  into  one  overall  loss  coefficient.  In  RTOD  the 
following  expression  from  reference  3  is  used  for  the  disk 
friction  loss  (eq.  (B 1 29)): 

0.02125  PiU\r\ 
df_  gcMotPMrj'n)0-2 

Rotor  tip  clearance  loss.  -  The  following  expression 
from  reference  3  is  used  in  RTOD  to  calculate  the  rotor 
tip  clearance  loss  (eq.  (B130): 


Z-ci  — 


r 4, tip) 

(1  —  ^.hub^^.iip) 


The  formulation  assumes  that  the  fractional  loss  due  to 
clearance  is  equal  to  the  ratio  of  clearance  to  passage 
height  at  the  rotor  exit. 


figure  3.  -  Geometry  input  variables. 


the  rotor)  such  that  the  RTOD  prediction  matches  two 
performance  parameters  (usually  flow  and  work  or  flow 
and  efficiency)  at  a  particular  pressure  ratio.  The 
parameters  to  be  matched  can  either  be  design-point 
values  obtained  from  a  design-point  program  or  test 
data. 


Program  Output 

The  program  output  is  in  either  SI  or  U.S.  customary 
units,  corresponding  to  the  chosen  input  units.  The 
output  consists  of  a  printout  of  the  title  card,  all  input 
data,  and  one  of  three  chosen  levels  of  output  as 
described  in  this  section.  The  RTOD  data  printout 
corresponding  to  the  input  data  of  figure  4  is  shown  in 
figure  5. 

Short  Output 

The  short  output  consists  of  a  tabulation  of  12  overall 
parameters  for  each  calculated  operating  point.  Figure  6 
shows  the  short  output  corresponding  to  the  input  data  of 
figure  4.  Note  that  the  rotor  chokes  and  that  the  value  of 
stator  exit  critical  velocity  ratio  ( V/  Fcr)|  at  rotor  choke  is 


printed  out.  Also,  the  rotor  limit  load  point  is  identified. 
The  parameters  in  the  short  output  are  as  follows  (with 
the  corresponding  equations  in  parentheses): 

(1)  Equivalent  flow  (B148) 

(2)  Equivalent  work  (B150) 

(3)  Actual  flow  (B44) 

(4)  Actual  work  (B 131) 

(5)  Equivalent  power  (B 151) 

(6)  Equivalent  torque  (B152) 

(7)  Equivalent  Po/ps  (B 153) 

(8)  Equivalent  p(>/pi  (B 154) 

(9)  Blade-jet  speed  ratio  (B162) 

(10)  Work  factor  (B161) 

(11)  Static  efficiency  (B134) 

(12)  Total  efficiency  (B 1 35) 

Medium  Output 

The  medium-output  option  gives  the  same  tabulation 
of  overall  parameters  as  the  short  output  and  also  gives 
tabulated  velocity,  temperature,  pressure,  and  flow  angle 
information  at  stations  0  to  5.  At  stations  3  to  5  the 
information  is  presented  for  both  the  absolute  and 
relative  frames  of  reference,  and  at  stations  4  and  5  (rotor 
exit)  the  information  is  presented  for  the  hub,  mid.  and 
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lip  sectors  only.  If  an  even  number  of  sectors  is  specified, 
the  midpoint  sector  is  taken  to  be  the  sector  above  the 
mean  radius.  The  medium-output  option  also  presents 
additional  information  for  stations  1,  3,  and  5.  For  the 
most  part  this  information  is  self-explanatory,  except  for 
these  items: 

Station  l 


stream  flow  (pO|  (equal  to  1.0  if 
p{/ p{)  was  specified;  less  than  1.0  if 
ei[)  was  specified). 

BL1  blockage,  defined  by  ratio  of  vane 

trailing-edge  tangential  thickness  to 
vane  trailing-edge  spacing 

Station  3 


EBR3D 


PIT/ POT 


V/VCR  AT 
MAX  FLOW 


three-dimensional  kinetic  energy  loss 
coefficient,  as  specified;  or,  if 
total  pressure  ratio  was  specified, 
calculated  from  equation  (B3) 
total  pressure  ratio,  p[/ pfr  (input 
parameter  PR  10),  as  specified;  or  if  ?3D 
was  specified,  calculated  from  equation 
(B4) 

critical  velocity  ratio  (V/V^ r)(  that 
produces  maximum  value  of  main- 


REL.  T.  PR.  relative  total  pressure  after  incidence 
AFTER  INCID.  loss,  a  fictitious  pressure  that  shows  the 
LOSS  effect  of  incidence 

U3  rotor  inlet  tip  speed 


Figure  7  shows  the  medium  output  for  a  (V/Va)s  of  0.6 
only  (for  the  example  input  of  fig.  4,  except  that 
OPT  =  I).  The  output  shown  is  produced  for  every  value 
of  ( V/  Fcr)(  in  its  specified  range. 


TABLE  I.  -  INPUT  VARIABLES 


— 

Input 

parameter 

Description  and  units 

_  .... 

Variable 

type3 

Default 

value 

Overall  input  variables 

TITLE 

Title  card  (up  to  80  alphanumeric  characters) 

- 

-- 

IUNTS 

Type  of  units  (0  for  U.S.  customary,  1  for  SI) 

I 

0 

TOT 

Total  temperature  at  station  0,  F  (°R) 

POT 

Total  pressure  at  station  0,  N/cm2  (psia) 

G 

Specific  heat  ratio 

R 

Gas  constant,  J/kg-K  (ft-lbf/lbm-°R) 

XMU 

Gas  viscosity,  kg/m-s  (ibm/ft-s) 

— 

OPT 

Output  option: 

0  -  short  output 

1  -  medium  output 

2  -  long  output 

0 

Stator  input  variables 

NSTV 

CLFR1 


Number  of  stator  vanes 

Stator  total  clearance  (hub  and  tip)  as  fraction  of  passage  height  at 
trail ing  edge 


0.0 


AL0 

AL1 

RO 

R 1 

HS0 

HS1 

THl 

I0PTS 


Flow  inlet  angle  at  station  0,  deg 
Stator  exit  blade  angle  at  station  1,  deg 
Radius  at  station  0,  cm  (in.) 

Radius  at  station  1,  cm  (in.) 

Passage  height  (including  clearance)  at  stator  inlet,  cm  (in.) 
Passage  height  (including  clearance)  at  stator  exit,  cm  (in.) 
Vane  trailing-edge  normal  thickness,  cm  (in.) 

Option  for  stator  loss  calculation: 

1  -  user  specifies  pj/pp 

2  -  user  specifies  633 


PR10 
E  BR  3D 
VRSTRT 

VRST0P 


Specified  ratio  of  Pj/Pq  (input  only  if  I0PTS  =  1) 
Specified  value  nf  ejp  (input  only  if  I  OPTS  =  2) 

Initial  value  of  (V/V.^)  station  1  (must  be  subsonic) 

Final  value  of  (V/V  )  (must  be  greater  than  VRSTRT) 

Cr  1 

Incremental  change  in  (V/V  )  between  VRSTRT  and  VRST0P 


R  0.98 

I  0.055 


DELVR 


TABLE  1.  -  Concluded. 


Input 

parameter 

Description  and  units 

Variable 

typea 

Default 

value 

Rotor  input  variables 

NRBl 

Number  of  full  rotor  blades  extending  from  inlet  to  output 

R 

I SPL 

Splitter  blades  (0  for  no  splitters,  1  for  splitters) 

I 

0 

RPM 

Rotative  speed,  rpm 

AlR3 

Rotor  inducer  sweep  angle  (0.0  for  radial  blade),  deg 

0.0 

R3 

R6WD 

R4HU8 

Radius  at  station  3,  cm  (in.) 

Rotor  "b"  width  at  station  3,  cm  (in.) 

Rotor  exit  hub  radius,  cm  (in.) 

R4TIP 

Rotor  exit  tip  radius,  cm  (in.) 

CL  4 

Clearance  height  at  rotor  exit,  cm  (in.) 

NSCT4 

Specified  number  of  equal  height  sectors  at  rotor  exit  (25  max) 

i 

0P8TA4 

Rotor  exit  blade  angle  distribution  option 

0  -  user  '..reifies  exit  blade  angle  for  each  sector 

1  -  program  linearly  interpolates  for  exit  blade  angle  at  each  sector 
from  input  BTAHU8  and  BTATIP 

I 

1 

8ETA4 

Rotor  exit  blade  angle  for  each  sector  (hub  to  tip;  input  only  if 

0PBTA4  =  0),  deg 

R 

BTAHU8 

Rotor  exit  blade  angle  at  hub  (input  only  if  0PBTA4  =  1),  deg 

R 

BTATIP 

Rotor  exit  blade  angle  at  tip  (input  only  if  0PBTA4  =1),  deg 

R 

0PTH4 

Rotor  exit  normal  thickness  distribution  option; 

0  -  user  specifies  normal  thickness  for  each  sector 

1  -  program  linearly  interpolates  for  normal  thickness  at  each  sector 
from  input  TH4HUB  and  TH4T1P 

i 

1 

TH4 

Rotor  exit  normal  thickness  for  each  sector  (hub  to  tip;  input 
only  if  0PTH4  s  0),  cm  (in.) 

1 

TH4HUB 

Rotor  exit  normal  thickness  at  hub  (input  only  if  0PTH4  =  1),  cm  (in.) 

TH4TIP 

Rotor  exit  normal  thickness  at  tip  (input  only  if  0PTH4  =  1),  cm  (in.) 

X.<R 

Rotor  viscous  loss  coefficient  for  each  sector  (hub  to  tip) 

PNEGI 

Power  in  incidence  loss  expression  (for  negative  incidence) 

2.50 

PPOSI 

Power  in  incidence  loss  expression  (for  positive  incidence) 

1.75 

P4RDPC 

Incremental  percentage  reduction  in  pressure  at  rotor  exit  after 
rotor  choke 

10.0 

dI  denotes  integer;  R  denotes  real. 
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8DATT  T0T=518.7.  F'0T  =  14.696.  G=l.4,  R=53.35>  XMU=1 .216E-5. 

NSTU  - 1 4  .  ,  AL0=45..  AL1=72.0.  R0=2.8B8.  Rl=2.345.  KS0=.577l>  HSU. 5771. 
IORTSM  »  PRIO*  .987 .  VRSTRT=.30.  URST0P=1.20.  DFLVR=.05.  THU. 02293. 

NRBL  =  1 1 . .  ISfL  =  l .  RF'N=29550 . .  R3=2.295»  RBHD=.577I.  R4HUB=.572>  R4TIP=1.646> 
C  L4=  .  0 1  45 »  NSCT4=11 .  0F'BTA4  =  1.  BTAHUB  =  -32 . 50 .  BTAT  IP=-56 . 70 .  0PTH4=1. 
TH4HUB=.07t,  TH4T  1P=  . 038 .  XhR  =  ll*.240,  0F’T=0  8 END 


Figure  A.  -  Example  1  data  input. 


Long  Output 


The  long  output  is  similar  to  the  medium  output, 
except  that  for  stations  4  and  5  the  tabulated  inputs  are 
presented  for  each  sector.  Also,  additional  rotor  loss 
information  is  presented.  For  the  most  part  this 
information  is  again  self-explanatory,  except  for  the 


following  items: 
EQUIVALENT  SPEED 
SPECIFIC  DIAM. 

SPECIFIC  SPEED 


defined  by  equation  (B149) 
specific  diameter:  defined  by 
equation  (B160) 
defined  by  equation  (B159) 


The  long  output  also  presents  a  further  tabulation  of  the 
following  items  for  each  sector  at  stations  4  and  5: 


Station  4 


R4 

BL4 

W/WCR  AT 


radius  associated  with  midpoint  of 
sector 

blockage  for  each  sector;  defined  by 
ratio  of  free  flow  area  at  station  4  to 
free  flow  area  at  station  5 
value  of  (W//ITcr)J  that  produces 
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IL)NTS=  0 

» 

T0T=0.5187£  03 

R  * 

P0T=0. 1470E  02 

RSIA 

G=0 . 1400E  01 

» 

XMU=0.1216E-04 

LBh/FUS* 

R=0 .5335E  02 
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Figure  5.  *  RTOD  data  printout  (corresponding  to  fig.  4  data  input). 
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Figure  6.  -  RTOD  short  output  for  figure  4  data  input. 
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Figure  7.  -  RTOO  mediun  output  at  tV'Vcr>  *  0.6  tor  figure  4  data  input. 


MAX  FLOW  maximum  (choking)  value  of  (p  WO4 
(always  less  than  1.0) 

U4  rotor  sector  speed 

Station  5 

AEROD.  ETA  sector  aerodynamic  static  efficiency; 

defined  by  equation  (BI27) 

AEROD.  ETAT  sector  aerodynamic  total  efficiency; 

defined  by  equation  (B128) 

WM(K)/WMTOT  fraction  of  total  weight  flow  passed 
by  each  sector 

Figure  8  shows  the  long  output  for  a  ( V/ Fa),  of  0.6  only 
(for  the  example  input  of  fig.  4,  except  that  OPT  =  2). 


The  long  output  shown  is  produced  for  every  value  of 
(V/Vcr)i  in  its  specified  range. 


Error  Messages 

All  DO-loops  in  RTOD  have  appropriate  error 
messages  that  are  printed  out  if  convergence  is  not 
achieved.  The  specified  convergence  criteria  are  defined 
by  the  ratio  of  the  absolute  value  of  the  new  value  minus 
the  old  value  to  the  old  value  and  are  usually  required  to 
be  less  than  1x10  4.  These  arc  light  tolerances,  but  no 
DO-loop  has  ever  failed  to  converge  for  physically 
possible  (low. 
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Figure  8.  *  RTOD  long  output  at  <V(  V 


*  0.6  for  figure  4  data  input. 


Main  Program 

The  user  can  specify  geometries  and  operating  condi¬ 
tions  that  have  no  flow  solutions.  For  low  flow  rates 
(small  specified  values  of  V/  Va  at  station  1 )  it  may  not  be 
possible  to  pass  the  flow  at  station  4  while  maintaining 
radial  equilibrium.  In  a  real  turbine  the  flow  would  be 
passed  but  with  separated  or  reversed  flow  in  portions  of 
the  rotor  exit.  However,  RTOD  cannot  account  for  that 
situation,  and  if  it  occurs,  the  following  error  message  is 
printed  out: 

STATION  4  HAS  NOT  CONVERGED  IN 
25  ITERATIONS  FOR  (V/VCR)1  =  x.xxxxx  ±xx 

ROTOR  WEIGHT  FLOW  =  xxxxxx  ±xx, 

STATOR  WEIGHT  FLOW  =x.xxxxx  ±xx 

After  this  message,  program  execution  starts  over  with 


the  next  specified  value  of  stator  exit  critical  velocity 
ratio. 

Subroutine  MAXVAL 

Subroutine  MAXVAL  calculates  the  maximum  value 
of  density  times  velocity  pV  and  the  associated  critical 
velocity  ratio  at  the  stator  and  rotor  exits  (stations  1  and 
4).  The  input  points  must  define  a  curve  that  shows  a 
maximum  somewhere  in  the  specified  interval  (other  than 
at  the  end  points).  If  unrealistically  large  values  of  stator 
or  rotor  loss  are  specified,  the  shape  of  the  generated 
curve  may  not  exhibit  the  required  characteristic.  If  this 
happens,  the  following  error  message  is  printed  out: 

WARNING -THE  CURVE  TO  BE  MAXIMIZED 
IN  SUBROUTINE  MAXVAL  IS  NOT  OF  THE 
PROPER  FORM 
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STATION  1 


PRIMARY  weigh:  flow  = 

0.6706E  00 

(LBM/S) 
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PRIMARY 
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=  Q.0000 

(LBN/S) 

ACTUAL 

PRIMARY 

FLOW  ANCLE 
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Figure  8.  -  Concluded. 


Program  execution  continues  after  this  message.  The 
value  of  pV  used  is  the  greater  of  the  two  at  the  interval 
end  points.  The  program  user  should  reexamine  and 
change  the  specified  stator  or  rotor  losses  if  this  error 
message  occurs. 


Example  Problems 

A  radial  turbine  that  was  tested  in  air  at  pressure  ratios 
higher  than  choking  is  analyzed  in  this  section,  and  the 
results  are  compared  with  experimental  data,  together 
with  the  results  from  the  program  of  reference  1.  Also, 
the  same  turbine  is  analyzed  for  various  values  of  rotor 
back-sweep  to  illustrate  the  potential  off-design 
performance  improvements  in  high-work-factor  turbines. 


Example  1  -  Comparison  with  Test  Data  and 
Reference  1  Code  Prediction 

Figure  9  shows  a  comparison  of  calculated  versus 
experimental  mass  flow  at  100-  and  60-percent  design 
speeds  for  the  11.66-cm  (4.59-in.)  diameter  turbine  of 
reference  5.  The  test  data  are  previously  unpublished  air 
data  that  were  used  because  they  covered  a  wider  range  of 
pressure  ratio  (beyond  rotor  choke)  than  the  published 
argon  data.  Both  the  RTOD  and  reference  1  code  predic¬ 
tions  were  matched  with  data  at  the  design  equivalent 
total-to-static  pressure  ratio  of  1.54.  The  RTOD  results 
more  closely  predict  the  experimental  mass  flow  rates  at 
high  pressure  ratios.  At  both  100-  and  60-percent  design 
speeds,  RTOD  predicts  an  approximately  1  percent 
greater  mass  flow  ratio  than  the  reference  1  code.  The 
improved  mass  flow  prediction  is  attributed  to  the  more 


Total  efficiency,  rf  Mass  ,low  rale  Percent  01  desi9n 


Figure  9.  -  Mass  flow  rate  comparison. 


accurate  geometry  and  How  modeling  at  the  rotor  trailing 
edge.  Both  the  RTOD  and  reference  1  code  total  effi¬ 
ciency  predictions  agree  closely  with  test  data  as  shown  in 
figure  10. 

Although  the  RTOD  mass  How  predictions  are  an  im¬ 
provement  over  the  reference  1  code  predictions,  they  still 
underpredict  choking  weight  flow  .  At  100-percent  design 
speed  the  RTOD  prediction  is  1.63  percent  below  test 
results;  while  at  60-percent  speed  it  is  2.28  percent  low  . 
This  discrepancy  could  be  caused  by  the  basic  rotor  loss 
modeling  scheme.  Rotor  viscous  loss  is  taken  to  be  a 
function  of  rotor  relative  inlet  and  exit  velocities  squared 
(eq.  (B87)),  multiplied  by  a  loss  coefficient  (input  param¬ 
eter  XKR)  that  is  assumed  to  remain  constant.  This  as¬ 
sumption  forces  the  calculated  maximum  mass  flow  rate 
per  unit  area  to  occur  at  low  values  of  rotor  exit  relative 
critical  velocity  ratio  (( H7/  W'cr)4  =  0.87  in  example  1). 
Proper  variation  of  the  rotor  loss  coefficient  with  relative 
flow  velocity  might  improve  the  choking  mass  flow 
prediction. 


Example  2  -  Effects  of  Rotor  Back-Sweep 

To  illustrate  the  effects  of  rotor  back-sweep,  the 
example  1  turbine  was  analyzed  with  assumed  back- 
sweep  angles  of  15°  and  30°.  The  predicted  total 
efficiencies  are  shown  in  figure  11.  For  high  values  of 
work  factor  (the  region  where  modern  radial  turbines  are 
likely  to  operate)  the  falloff  in  efficiency  is  less  with  back- 
sweep  than  for  a  radial-bladed  rotor  (because  the 
incidence  is  lower). 


National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio,  March  25,  1983. 
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Appendix  A 
Symbols 


A  area,  m-  (ft-) 

BL.4  ratio  of  blocked  area  at  station  4  to  free  flow 
area  at  station  5 

CLFR1  ratio  of  total  clearance  area  to  total  passage 
area  at  station  1 

Df,  hydraulic  diameter,  m  (ft);  defined  by 

eq.  (B46) 

Ds  specific  diameter;  defined  by  eq.  (B160) 

eio  three-dimensional  kinetic  energy  loss 

coefficient;  defined  by  eq.  (B3) 

/  friction  factor 

gc  force-mass  conversion  constant,  1  (32.174 

(lbm-ft  )/(lbf-s2)) 

h  height,  m  (ft);  or  enthalpy,  J/kg  (Btu/lbm) 

J  dimensional  constant,  1  (778  ft-lbf/Btu) 

A'  rotor  viscous  loss  coefficient 

k  sector  number 

L  kinetic  energy  loss,  J/kg  (Btu/lbm) 

/  length,  m  (ft) 

M  Mach  number 

N  number  of  sectors  at  station  4 

Nbi  number  of  rotor  blades  (including  splitters) 

Ns  specific  speed;  defined  by  eq.  (B159) 

n  incidence  exponent 

P  power,  W  (hp) 

p  pressure,  N/m2  (psia) 

Q  parameter  defined  by  eq.  (B166) 

R  gas  constant,  J/kg-K  (ft-lbf/!bm-°R) 

Re  Reynolds  number;  defined  by  eq.  (B47) 

r  radius,  m  (ft) 

T  temperature,  K  (°R) 

TRQ  torque,  N-m  (ft-lbf) 

U  tangential  velocity,  m/s  (ft/s) 

V  velocity  in  absolute  reference  frame,  m/s 

(ft/s) 

W  velocity  in  relative  reference  frame,  m/s  (ft/s) 

WRK  work,  J/kg  (Btu/lbm) 

w  flow  rate,  kg/s  (Ibm/s) 

a  fluid  flow  angle  in  absolute  frame  of  refer¬ 

ence;  or  blading  angle,  deg 

d  fluid  flow  angle  in  relative  frame  of  reference; 

or  blading  angle,  deg 
y  ratio  of  specific  heats 

6  parameter  defined  by  eq.  (BI56) 


<  parameter  defined  by  eq.  (B157) 

f|jn,  parameter  defined  by  eq.  (B73) 

i)  efficiency 

0cr  parameter  defined  by  eq.  (B155) 

n  viscosity,  kg/m-s  (Ibm/ft-s) 

v  loading  coefficient  (blade-jet  speed  ratio); 

defined  by  eq.  (B162) 
p  density,  kg/m2  (lbm/ft2) 

work  factor;  defined  by  eq.  ( B 1 6 1 ) 

(2  rotational  speed,  rpm 

Subscripts: 

act  actual 

aed  aerodynamic 

av  average 

bl  blading 

ck  choked 

cl  clearance 

cr  critical 

df  disk  fraction 

equ  equivalent 

/  friction 

hub  hub 

i  incidence 

id  ideal 

mx  maximum 

opt  optimum 

p  primary 

r  radial  component 

ref  reference 

shr  shroud 

std  standard 

tip  tip 

tot  total 

u  tangential  component 

x  axial  component 

0,1,2,  station  numbers  defined  by  fig.  1 

3,4,5 

Superscripts: 

(  )'  total  conditions  in  absolute  reference  frame 

(  )' '  total  conditions  in  relative  reference  frame 

(  )  mass  averaged 


Appendix  B 
Equations 


Main  Program  (MAINP) 

The  analysis  of  the  flow  through  the  turbine  starts  at 
the  stator  exit  (station  1)  for  each  specified  value  of 
(V  l  Fcr)  .  Flow  conditions  are  then  calculated  at  station  0 
and  at  stations  2  to  5.  The  procedure  is  described  in  this 
appendix. 

Station  1  —primary  flow.  -Total  temperature  at 
station  1  is  assumed  to  be  equal  to  that  at  station  0.  Thus 

T[,P=n,p  (Bn 


Both  methods  of  calculating  the  total  pressure  drop  are 
equivalent  to  specifying  an  overall  viscous  loss  that  takes 
into  account  the  vane  profile,  endwall,  and  secondary 
flow  losses. 

From  the  isentropic  relationships 


The  maximum  value  of  primary  flow  per  unit  area 
[(pFJj  ]  and  the  value  of  (V/Vcr),  where  that 
maximum  occurs  l(V/Vct)l  are  calculated  by 
subroutine  MAXVAL.  If  the  total  pressure  drop  through 
the  stator  (station  0  to  station  1)  is  modeled  by  specifying 
the  total  pressure  ratio  p{/po,  (pV) ,  x  occurs  at 
(V/VCT)[  mx  =  1.0.  If  the  total  pressure  drop  is  modeled 
by  specifying  ejp-  it  occurs  at  (k7F'cr)  x  <  1.0.  If 
p\/pb  is  specified,  total  pressure  at  station  1  is  calculated 
from 


(B6) 


(B7) 


and  the  three-dimensional  kinetic  energy  loss  coefficient  For  values  of  ( V/  Vcr)i  less  than  ( V/  mx  the  primary 
?3D  is  obtained  from  flow  exits  at  the  specified  stator  exit  blade  angle  oq.  Thus 


'v\,p  =  P\,pv\,pA\tt  -CLFR1)  cos  oq>p  (B9) 

For  values  of  ( V/  Kcr)]  greater  than  ( V/  ^r),,  (choked 
flow)  the  exit  angle  is  decreased  to  maintain  maximum 
primary  flow  rate.  The  choked-flow  exit  angle  is 


ul,p,ck  —  COS 


L.^>„ 


pA\(\ -CLFR1) 


If  ejD  is  specified, 


(B4) 


The  radial  and  axial  components  of  the  primary  flow  are 
then  defined  by 

^u,\,p=  ^i.p  sin  ( B 1 1 ) 

vr,\,p=  h.pcos  al<p  (B12) 

Mach  number  is  calculated  from 


u 


r 


If  the  radial  component  of  Mach  number  cos  ct\  p) 
is  greater  than  1 .0,  the  specified  ( V/  l/cr)1  is  reduced  until 
M\<p  cos  ot  \  p  =  1.0  (defined  as  limit  loadfing). 

Station  1— clearance  flow.  -  The  clearance  flow  is 
assumed  to  have  the  same  total  and  static  conditions  of 
temperature  and  pressure  as  the  primary  flow: 


Since  aij  is  uniquely  defined  by  equation  (B24),  it 
cannot  be  changed  to  maintain  the  maximum  clearance 
weight  flow  at  values  of  ( V  / 1  cr),  greater  than 
( VI  ^-r),  |  itu.  Clearance  weight  How  therefore  decreases 
above  (K/Pcr)|  (  iiix. 

From  here  on,  all  conditions  at  station  I  that  are  the 
same  for  both  the  primary  and  clearance  flows  ( T,  T' ,  p, 
p' ,  V,  Fcr,  and  p)  will  have  the  subscripts  p  or  cl  deleted 
and  will  be  identified  only  by  the  subscript  1. 

Station  0.  -  Conditions  at  station  0  are  evaluated  from 
the  specified  geometry,  inlet  angle,  and  gas  total 
conditions  and  by  using  mass  continuity  between  stations 
0  and  1 .  The  ratio  ( V/  Fcr)(|  is  determined  iteratively  from 


(M’i,p+  w,t,c|)/?7o 


7-i.ci  =  r[p 

(B 14) 

(  V\  - 

V^r/O 

r  -v  - 1  / 

i/Vl'B  i) 

7"i  ,ct =  T\  p 

(B15) 

A0  COS  aoPoFo.crl  ■  7  +  l^ 

Pi, cl  ~P\.p 

(B16) 

(B26) 

Pi, cl  ~P\,p 

(B17) 

Also 

This  implies  that 

f  y  \  2~|  >/(T  “  ■> 

,^/oJ 

^,|,cl.cr=  ^l.p.er 

(B18) 

*>-*[■ -*tt( 

(B27) 

h,d=h,p 

(B19) 

/  V  \  /  V  \ 

ro-r6[i-;+;i 

taj 

(B28) 

V^A.c=  w/i.p 

(B20) 

(  V  \  (  V  \ 

p°=  £t0 

(B29) 

(B21) 

l,p,mx 

P|.Cl=P|,p 

(B22) 

v  _(  V\  (2yRgJ»\'> 

(B30) 

The  flow  angle  «|tC| 

and  the  flow  rate  W|  c(  are 

yu,0  =  >o  sin  a0 

(B31) 

determined  iteratively  from  conservation  of  moment  of 

tangential  momentum 

between  stations  0  and  1.  The 

Vr.O  =  V0  COS  a0 

(B32) 

tangential  component  of  velocity  is  obtained  from 

Station  2.  -  The 

total  temperature 

at  station  2  is 

,/  '0^.0 
»Vi.d- 

(B23) 

assumed  to  equal  that  at  station  1 

'  1 

72-77 

(B33) 

The  clearance  flow  angle  and  flow  rate  thus  are 


al,cl  =  sin  ~  1 


(B24) 


and  tangential  momentum  is  conserved  between  stations 
1  and  2.  The  tangential  velocity  is  a  mass-averaged  value 
of  the  primary  and  clearance  flows  at  station  1: 

V\ 

y„  i  =  (  H'|  sin  r*|  „+  W|  c|  sin  <V|  c|)  (B34) 

w’l./>  +  w’l.d 


"l.cl-Pl.d ^l.d^l  CLFR1  cos  a,td 


(B25) 


For  no  clearance  flow  this  reduces  to 

vu.2=  v\  sin  «i,p=  ^u.i  (B35) 

By  using  an  initial  guess  for  Vr  i2,  the  conditions  at  station 
2  are  evaluated  iteratively  as  follows: 


Re2  = 


P2  V2Dh 
P 


/2-3  =  SlZ-Ll— 

COS  a2_3 


(B47) 


(B48) 


y2=(V2r.2+V2u/ 


V2 


(B36)  Friction  factor  /is  defined  as  follows: 

/  =0 . 2  4  for  Re2 <  100  (B49) 

(B37) 

/=  for  100  <Re2<  3000  (B50) 


P2=Pi 


P2=P2 


7+lV^cr/2_ 


•>''(7-0 


a2  =  COS  - 


W2  =  P2V2^2  cos  a2 


(B38) 

(B39) 

(B40) 

(B41) 

(B42) 

(B43) 


The  assumed  Vr2  is  updated  until  continuity  is  satisfied: 
that  is. 


/=  3.3368  x  10-7  ReJ  2596  for  3000  <Re2<  3700  (B51) 

7=0.0014+  for  3700  <Re2  (B52) 

Re2 

The  change  in  total  pressure  and  the  total  pressure  at 
station  3  are  obtained  from 


y ,  Vhinvj 
^  2  Dm 

(B53) 

pi=p2-Ap’ 

(B54) 

Assuming  that 

K1 

11 

KT 

(B55) 

r2  =  r\ 

(B56) 

and 

M'2=lvl,p+vvl,cl  (B44) 

Station  3.  -  The  total  pressure  drop  between  stations  2 
and  3  is  obtained  from  the  laminar  and  turbulent  equa¬ 
tion  for  flow  between  parallel  plates  by  using  an  average 
flow  angle  between  stations  2  and  3.  Furthermore, 
moment  of  tangential  momentum  is  conserved 
(r2  yu,2  =  ^3  Vu,})-  From  an  initial  guess  of  03  an  average 
turning  angle,  a  hydraulic  diameter,  a  Reynolds  number, 
and  a  flow  path  length  are  calculated: 


on  +  c*3 

«2-3=-y- 


(B45) 


^.3  = 


r2/>3 


(B57) 


all  other  conditions  at  station  3  are  evaluated  from 


^3 


sin  c*3 


(B58) 


,  /  7\  \7/(7  -  I) 

P3=/?3^~j  (B6()) 


Dh  =  2h\  cos  a2-3 


(B46) 


IS 


,,  f /r,2PlAl 

'  r,3  —  , 

P3^3 

(B62) 

Finally,  a  new  exit  angle 

a,"an"(fe) 

(B63) 

is  calculated  and  compared  with  the  previously  used 
value.  This  procedure  is  repeated  until  successive  values 
of  a3  are  within  a  specified  tolerance. 

Relative  conditions  at  station  3  are  calculated  from  the 
following  equations: 

W  u.3  =  f  «,3  ~  t/j 

(B64) 

MV.3=  Vr.\ 

(B65) 

=  (  W2  ,4-  )l/2 

-  u.  3  r,  3' 

(B66) 

(B67) 

3a 

II 
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UJ  - 

/HI  *H 

»  1  5 

1 

(B68) 

(B69) 

(B70) 

The  rotor  incidence  angle  (from  optimum)  is  calculated 
as  follows:  the  optimum  Vu i3  is  obtained  from  the 
Wiesner  slip  factor  correlation  (ref.  4),  which,  when 
adapted  to  a  turbine  geometry,  takes  the  form 


..  ^4, hub  +  ^4, up 

f4.ai-  2 - 

(B72) 

1 

f''m  ^(8.16  cos  <«b|/.Vhlj 

(B73) 

The  last  term  in  the  numerator  of  equation  (B71)  is  a 
correction  factor  that  is  applied  only  if  the  ratio  r^/ry  is 
greater  than  f|im.  Then 


^  u,3*opt  ""  ^w,3,opi  ^3 

(B74) 

(B75) 

«3.(  =  ^3  -03.opt 

(B76) 

The  loss  in  total  pressure  due  to  incidence  is  calculated 
as  follows:  the  incidence  loss  is  modeled  as  in  reference  2 
by 

-cos*a3 ,  ) 

L,~  2  g(J 

(B77) 

where  the  exponent  n  is  given  by  the  input  parameters 
PNEGI  and  PPOSI  for  negative  and  positive  incidence, 
respectively.  The  loss  due  to  incidence  can  also  be 
expressed  as  a  difference  between  ideal  and  actual 
velocity: 

W]-  W\.=  W]{\  —  cos"a3t/) 

(B78) 

Solving  for  W2  .  gives 

^3  i  =  ^3  cos”a3.i 

(B79) 

The  relationship 


^u,3,opt  ~ 


( 


1  -  Vcos  ab|/(/Vb|)0.7][i  -[(r4av/r3-£|jm)/(l  -eiim)]3; 


1  -  tan  aw/tan  a3 


) 


(B71 ) 


implies  that 

Pi 

P  3,/  = - - 

(y- \)W*  ‘1>/(T,-I) 

1_  2yglJtT/ 


(B80) 


(B81) 


Substituting  equation  (B79)  for  H'j  j  then  gives 


station  4  p'i  is  obtained  from  p 4^  by  considering  the 
combined  effects  of  incidence  and  rotor  viscous  loss. 

The  value  of  maximum  flow  rate  per  unit  area 
( p  ff)4,mx,  as  well  as  the  value  of  ( W/  Wa)t  where  it  occurs 
\(W/  mJ,  is  calculated  by  subroutine  MAXVAL. 
The  total  relative  pressure  at  station  4  is  formulated  as 
follows:  the  incidence  loss  is  modeled  by  (see  previous 
section) 


^-(l  -cosna3-() 

2^7 


(B86) 


- _ -  (B82) 

(y-\)W\  cos«a3,'"|>/l>-1» 

2y  gcRTS 


Station  4.  -  Flow  conditions  at  station  4  are  calculated 
by  considering  the  effects  of  radius  change,  incidence, 
and  viscous  loss.  By  assuming  an  initial  value  of  static 
pressure  at  the  hub  sector  the  pressure  at  each  successive 
sector  is  obtained  from  simple  radial  equilibrium.  The 
total  mass  flow  of  all  sectors  is  then  compared  with  the 
mass  flow  at  station  1,  and  the  next  guess  of  hub-sector 
static  pressure  is  changed  accordingly.  This  is  repeated 
until  mass  flow  has  converged  within  a  given  tolerance. 
The  procedure  is  given  here  (with  calculations  performed 
for  each  sector  at  the  midpoint  radius). 

The  changes  in  total  temperature,  critical  velocity,  and 
ideal  relative  total  pressure  due  to  radius  change  are 
calculated  from 

TS=T!+m-Jui-u9  (B83) 


and  the  loss  due  to  viscous  friction  in  the  rotor  (rotor 
blade  profile,  end  wall,  and  secondary  flow  losses)  is 
modeled  as  in  reference  1  by 


Lf= 


K(W]  cos2a3,+  H^2) 

2g<7 


(B87) 


where  K  is  the  specified  rotor  viscous  loss  coefficient 
(input  parameter  XKR).  The  sum  of  these  two  losses  is 
expressed  as  the  difference  between  the  squares  of  an 
ideal  and  an  actual  velocity  to  give 


+  K(  W 2  cos2a3  ,  +  W (B88) 


This  equation  can  be  solved  for  WA 


^(l  -  cos ”a3,/)  -K tVj  cos2a3  ,  "l  1/2 


1  +K 


(B89) 


/  77\l/2 

^cr.4=^cr,3(^j  (B84) 

/  1) 

PiM=Pi\Y.)  (B85) 

Note  that  the  ideal  relative  total  pressure  at  station  4  is 
based  on  the  relative  total  pressure  at  station  3  and  not  on 
the  relative  total  pressure  due  to  incidence  loss  pjj.  At 
station  3,  pj j  was  calculated  for  information  purposes 
only.  As  will  be  shown,  the  relative  total  pressure  at 


with  fF2  obtained  from  the  relationship 


4, id 


^?r.4(7+  f,  (  Pi 

[  WJ 


(B90) 


For  a  given  value  of  static  pressure  the  relative  velocity  is 
thus  calculated  by  equation  (B89).  Also, 


74=74'- 


2yRgc 


(B91) 
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p4  = 


P± 

rt4 


(>  - 1) 


(B92) 


(B93) 


The  static  pressure  at  each  radial  sector  above  the  hub 
sector  is  obtained  from 


dp  _  p4jk-  1)|[F(*~  1)„.4+  m)M<4]/2)2 

dr  gcr(k  - 1)4 

p4(k)  =p4(k-  l)+  ^Ar 


(B94) 


(B95) 


It  is  not  allowed  to  exceed  1.0.  If  it  does,  the  value  of 
hub-sector  static  pressure  is  adjusted  until  Mx4  equals 
1.0.  That  value  of  hub-sector  static  pressure  is  also 
printed  out. 

The  absolute  conditions  are  calculated  from  the 
relative  conditions  by 


^.4  =  WU,4+t/4 

Kv,4  =  M\,4 

F4=(F2  +  y2  ) 1  2 

*  u, 4  v,4 

“4=,an'(^) 


(BI01 ) 
(B102) 

(B103) 

(B104) 


where  k  designates  the  given  sector  and  (k  -  1)  designates 
the  sector  immediately  below  it. 

For  ( W/  FF'cr)4  less  than  (W/Wcr)4mx  the  flow  follows 
the  specified  exit  angle  /34;  for  (M//M/cr)4  greater  than 
( W/  B/Cr)4  the  exit  angle  is  decreased  to  maintain  the 
maximum  flow  rate.  The  reduced  exit  angle  is  calculated 
from 


3  =  cos  -  1  T  1 

04  L  P4  W4A 5(  1  -  BL4)  J 


(B96) 


where  BL4  is  the  ratio  of  the  blocked  area  at  station  4  to 
the  free  flow  area  at  station  5.  Weight  flow  and  velocity 
components  are  calculated  from 


w4  =  p4WVl5(l  -BL4)  cos  04 
Wu4  =  W4  sin  04 
Wx  4  =  W4  cos  04 


(B97) 


(B98) 

(B99) 


If  the  rotor  is  fully  choked  (all  sectors),  the  program 
solves  for  that  value  of  stator  where  rotor  choke 

occurred.  Also,  the  maximum  value  of  station  4  hub- 
sector  static  pressure  for  rotor  choke  is  calculated  and 
printed  out.  Once  the  rotor  has  fully  choked,  the  hub- 
sector  static  pressure  is  incrementally  decreased  by  a 
specified  amount  (input  parameter  P4RDPC)  until  limit 
loading  is  achieved  (axial  Mach  number  equal  to  1.0). 
The  axial  component  of  Mach  number  is  calculated  from 


MxA  = 


J 

4t)< 

f  W  \2 

v  ^cr  /  4 

1  - 

1/2 


COS  04 


(B100) 


n-n* 


2  yRgc 


P4- 


P4 


(74/74) 


( B 1 05 ) 


(B106) 


Station  5.  -  At  station  5  the  flow  area  is  divided  into 
the  same  sectors  as  at  station  4.  All  calculations  shown 
below  are  performed  for  each  sector  at  the  midpoint 
radius.  The  total  temperature  and  tangential  momentum 
at  station  5  are  assumed  to  be  equal  to  those  at  station  4: 


^,5=^,4 

Equation  (B107)  implies  that 

^cr,5  =  ^cr,4 


(B107) 

(B108) 

(B109) 


From  an  initial  guess  of  05  the  flow  conditions  are 
established  iteratively  as  follows: 


sin  05 


\KrJi  Kr.S 


7-5=  7"5- (7-  1) 


J 1 

2yRgc 


,  ,  Ft  cos  05 \2 /O' •5P4^ '4  \ 

{~nr) 


(B 1 10) 


(Bill) 


(B 1 12) 


(B 1 13) 
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-['-(&)(ic)3 


y'lh  -  1) 


„  _  £i_ 

P5  /?r5 


_  VxAP4A5(l  -BL4) 
P5»5 


(B 114) 


(B 115) 


(B 116) 


A  new  exit  angle  is  then  calculated  and  compared  with  the 
previous  value. 


«5  =  tan  - 1 


(fe) 


(B 117) 


The  procedure  is  repeated  until  successive  values  of  a 5  are 
within  a  specified  tolerance. 

Performance  parameters.  -  The  performance  param¬ 
eters  are  calculated  as  follows:  for  each  sector  the  static 
and  total  ideal  enthalpy  change  from  station  0  to  station  5 
is  calculated  from 


(7- DJ 


(B 118) 


(B 1 1 9) 


Mass-averaged  values  of  ideal  enthalpy  change,  power, 
and  torque  are  established  from 


‘*’101=  E  wik) 

k  =  I 


(B120) 


Ahn<  = 


E  ±h(k)0.5wik) 

k  =  I 


E  Ah'  (Ar)0_,w(*) 

*  =  1 


(B121 ) 


(B122) 


wk„7  A/to-5 


TRQ  = 


wtoty  A/j6_5(60) 


(B124) 


The  change  in  total  enthalpy  across  the  rotor  (station  3  to 
station  5)  for  each  sector  is  calculated  from 


....  u,v„j-u,v„A 

— jfe — 


and  its  mass-averaged  value,  from 


(B125) 


A/jj-5  = 


2]  Ah'  (k)i.5w(k) 
*=! 


(B126) 


The  aerodynamic  efficiencies  for  each  sector  are 
calculated  from 


7aed  — 


(B127) 


(B128) 


Disk  friction  and  exit  tip  clearance  loss  per  unit  mass  flow 
are  obtained  from  (ref.  3) 


0.021 25P)U]r] 


(B129) 


(B130) 


Actual  total  enthalpy,  actual  power,  and  actual  torque 
across  the  rotor  are 


A/taci  =  hi-$~  Ldt~  Ls\ 


_  wto\J  A/tacl 
aa  550 


(B 1 3 1 ) 


(B132) 


60H-ri,,y  AhL. 


and  actual  overall  efficiencies  are 


Rslli  =  287.05  J/(kg-K)  (53.35  ft-lbf/lbm-°R)  (BI46) 


7std  =  1-40 


(BI47) 


(B134) 


(B135) 


Mass-averaged  static  and  total  temperatures  are 
calculated  by 


£  T(k)<w{k) 

k  =  I 


(B136) 


The  equivalent  parameters  (flow,  speed,  work,  power, 
torque,  total-to-static  pressure  ratio,  and  total-to-total 
pressure  ratio)  are 


ct|u~  s 


^equ  ~ 


WRK-=  if 


(B148) 


(B149) 


(B150) 


£  T'{k)<w(k) 
k  =  1 


(B137) 


p  _  ^actf 

'equ  -  - 

8^ 


(B 1 5 1 ) 


Static  pressures  at  the  hub  and  shroud  and  average  static 
and  total  pressures  are  obtained  from 


(d£\  = 
V  dr  /hub 


P(1)5^2(DM.5 
gA  1)4 


(dp\  p(N)iVHN)Uti 
\dr  Ahr  gfW)A 


^.hub-P'D,- (*)lubf 


P 5.av  — 


PS.hub  +  /*5.shr 


Pl.av  ~  Ps,a\  ^  ~  ^ 


V<7  -  I) 


(B138) 


(B139) 


(B140) 


(B 141) 


(B142) 


(B143) 


Equivalent  parameters  are  calculated  on  the  basis  of  the 
following  standard  conditions: 


7",',d  =  288.2  K  (5 18.7  °R) 

Pq,\  =  101  325  N/m2  (21 16.22  lbf/ft2) 


(B144) 


TRO  — 

*  *'Vequ  —  g 


(B152) 


/Po\  _  Tj  _  ( Tstd  ~  1  )J  A/tp-5  1  T«d/('K(d  I) 

\P$ ) equ  L  ^cr'Vstd^std^std  J 


(po\  r1  (Ts.d- 1)-/  a/iq- 

\PS/e  qu  L  ^cr7std^std^"std 


( 7sld  _  1  )J  M6-S  1  “  Wit-id  -  I ) 


(B153) 


where 


(^l)Rs,d7’S',< 


6=^ 


2  \7sld/(7sld  -  •) 


Mil 

<£) 


y/(y  -  1) 


(B154) 


(B155) 


(B156) 


(BI57) 


(B145) 


The  overall  dimensionless  parameters  of  specific  speed, 
specific  diameter,  work  factor,  and  loading  coefficient 
(blade-jet  speed  ratio)  are  defined  as  follows: 


r.  =  7T 


HsO(«] 


(B165) 


Ps.av 


Ps.av 
R  ^5,av 


<bis81  e=G^)(H 


(B166) 


/2jrfi 
\  60 


)(£) 


1/2 


aao-s)°-75 


2r-\(gcJ  A/iq-5)()  2? 
(‘iW/P5.av),/2 


If  the  total  pressure  at  station  1  is  calculated  by  specifying 
IOPTS=  1  (specifying  the  pressure  ratio  p\/po), 

(B159) 


Pl  = 


(B167) 


(B160) 


If  the  total  pressure  at  station  1  is  calculated  by  specifying 
lOPTS  =  2  (specifying  e-}o)< 


gtJ  A//j_5 
"1 


(B 1 61 ) 


Pi=A)[ 


I  -e3D-g 
(l-<?3D)(l-0 


] 


(B168) 


Uj 

(2glJ  A/x0-5>i/2 


(B162) 


Subroutine  MAXVAL 

Subroutine  MAXVAL  finds  the  maximum  value  of  a 
function  over  an  interval.  The  chosen  function  must  be  of 
a  form  such  that  a  maximum  exists  within  the  interval 
(other  than  at  the  end  points).  MAXVAL  is  used  to 
calculate  the  maximum  values  of  flow  rate  per  unit  area 
at  the  stator  exit  (pV)  |  and  rotor  exit  (pHO.*  sectors.  The 
actual  mass  flow  calculations  are  performed  by  the  called 
subroutines  ROVST  and  ROVROT  for  the  stator  and 
rotor,  respectively.  Given  the  interval,  the  function  to  be 
maximized,  and  an  initial  number  of  divisions  of  the 
interval,  MAXVAL  determines  the  function’s  maximum 
value  and  its  location  in  the  interval. 


Subroutine  ROVST 


Finally, 


7/(>  -  1) 


P\  =P\  J^l  -  1 

(BOK)!] 

(B169) 

Pl  RT\ 

(B 1 70) 

pK=p,  F, 

(B I 7 1 ) 

Subroutine  ROVROT 

Subroutine  ROVROT  calculates  the  value  of  density 
times  velocity  at  the  rotor  exit  sectors  (station  4)  for  the 
assumed  rotor  loss  model  and  a  specified  ( W/  (Tcr)4.  The 
equations  used  are  as  follows: 


Subroutine  ROVST  calculates  the  value  of  density 
times  velocity  at  the  stator  exit  (station  1)  for  the  assumed 
stator  loss  model  and  a  specified  (V/V„)r  The  formula¬ 
tion  is  as  follows: 


(B 1 63) 


(B164) 


( B 1 72) 


Tt=TS- 


(7-1 


(B173) 


The  static  pressure  is  calculated  from  equation  (B88), 
which  expresses  the  sum  of  the  two  rotor  losses  (incidence 
and  viscous  friction)  as  the  difference  between  the 
squares  of  an  ideal  and  an  actual  velocity 
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w-  -  W'2  =  B'2(  1  - 

4, id  4  3'* 

COS^Q^j) 

i  ncn 

+  K{  W\  cos-aj  |  +  W'2) 

(B174) 

PA=Pa\  1 

By  using 

L 

“  (^j)”  ‘ ,V>] 

(B 1 75 ) 

«te’ 

II 

~r 

Q. 

Hi-ri 

(B176) 

pV=piWi 

equation  (B174)  can  be  rewritten  as 

(izl 
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